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ABSTRACT 


The concept of a hypersonic research airplane is one which has 
been studied for several years. The project reported on in this 
paper involved the testing of a hypersonic beaded skin panel to 
failure. The primary interest was focused upon the buckling char- 
acteristics of the panel under pure compression with boundary con- 
ditions similar to those found in a wing mounted condition. Three 
primary phases of analysis are included in this report. These phases 
include: Experimental testing of the panel to failure; Finite element 

i 

structural analysis of the beaded panel with the computer program 
Nastran; A summary of the semiclassical buckling equations for the 
beaded panel under purely compressive loads. A comparison of each 
of the analysis methods is also included. 
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Cross sectional area of section defined by half pitch 
dimension, a, (see fig. 28) at the center cross section 
of the panel, inches 2 (cm-) 
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Bending or twisting rigidity, LBf*in (N-m) 
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Critical „axial compressive stress for bead instability, 
psi (N/nr) 
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Critical bendinq stress for bead instability, psi (N/m ) 

Critical axial buckling load for flat instability, 
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Moment of inertia per unit length, XY plane 
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Critical axial buckling load, LBf (N) 

Critical axial buckling load for general instability, LBf 
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Critical axial buckling load for diagonal instability, LBf 
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SUMMARY 


The concept of a hypersonic research airplane is, one which has 
been studied for several years. Early studies included research into 
new structural concepts with emphasis placed upon developing the best 
cost/weight efficiency, performance and reliability obtainable. 

As a part of NASA's continuing research into hypersonics, Dryden 
Flight Research Center has been laboratory testing an 85 square foot 

p 

(7.9m ) hypersonic wing test section of a proposed hypersonic research 
airplane. In tests performed to date on the wing test section, the 
structure has exceeded all expectations of strength and durability 

The project reported on in this paper has carried the hypersonic 
wing test structure project one step further by testing a single 
beaded panel to failure. The primary interest was focused upon the 
buckling characteristics of the panel under pure compression with 
boundary conditions similar to those found in a wing mounted condition. 
Three primary phases of analysis are included m the report. These 
phases include: Experimental testing of the beaded panel to failure; 
finite element structural analysis of the beaded panel with the com- 
puter program Nastran; a summary of the semiclassical buckling equations 
for the beaded panel under purely compressive loads. Comparisons be- 
tween each of the analysis methods is also included. 
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INTRODUCTION 


The Hypersonic Research Airplane 

The concept of a hypersonic research airplane (HRA) is one which 
has been studied for several years (ref. 1-10). Early studies included 
not only basic conceptual design, but also research into new structural 
concepts. This research has provided a portion of the technological 
base necessary for future hypersonic developments. 

Much of the research done has been devoted to theoretical analysis 
of various structural concepts which meet the requirements of a hyper- 
sonic airplane. Emphasis has been placed upon developing the best 
cost/weight efficiency, performance and reliability obtainable. Weight 
efficiency in high performance aircraft is a critical factor explain- 
ing the need for a weight efficient structure. 

One HRA concept studied by NASA is shown in figure 1. This pro- 
posed vehicle would cruise at Mach 8 for five minutes. It is a single 
place design with a wing span of 38 feet (11.58m), a length of 101 feet 

c 

(30.78m) and an estimated weight of 75,600 pounds (3.36x10 N). The 
wings and tail are hot radiating structures fabricated from super 
alloys. 


2 




Figure 1. Hypersonic research airplane configuration concept. 
Note • Dimensions are in feet and (meters). 
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The Hypersonic Wing Test Structure 


As part of NASA's continuing research into hypersomcs, Dryden 
Flight Research Center has been laboratory testing an 85 square foot 
(7.9m ) hypersonic wing test section, shown in figure 2, of the pro- 
posed HRA vehicle. The objectives of this program are to verify ana- 
lytical predictions, construction techniques, assembly techniques and 
in general to improve flight loads measurement technology. 

The hypersonic wing test structure (HWTS), shown in figure 3, is 
made from Ren£ 41 (with the exception of the lower leading edge heat 
shield panels which are TD Ni Cr) and is capable of operating with 
surface temperatures in excess of 1800° F (1250°K). The HWTS employs 
corrugated spar and rib webs and beaded skin panels. Aerodynamic 
smoothness is accomplished by attaching heat shields over the beaded 
panels. 

The^HWTS carries loads somewhat differently than do conventional 
aircraft. Bending loads normally carried by spars in conventional 
wing structures are instead carried by the beaded skin panels in the 
HWTS. Shear and torque are carried in much the same manner as in con- 
ventional wings. 




Figure 2 . Hypersonic wing test section of the proposed hypersonic research airplane . 
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Figure 3. The hypersonic wing test structure ( HWTS ) 
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transducers (DT's). The second technique used was the f^oird fringe 
(or grid shadow) technique (ref 12). The Moird fringe technique oro- 
vided deflection data for the entire panel versus the descrete measure- 
ment of the DT. The DT data also provided a check of the Moird data. 

A brief description of the Moird fringe technique is included in appen- 
dix C. 

The second phase of the project was the finite element structural 
analysis of the panel. The stresses, deflections, and buckling char- 
acteristics were calculated with Nastran (ref 13) a finite element 
structural analysis computer program. 

The third phase of the project included a summary of the semi- 
classical analysis previously done for the beaded panels. Semiclassi- 
cal equations of buckling strength were developed for beaded panels of 
the type used for this test in references 2, 3, and 4. 
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THE BEADED PANEL TEST SPECIMEN 


Beaded Panel Description 

The beaded panel concept meets the requirements of high strength 
and weight efficiency required for a hypersonic airplane. The panel, 
as shown in figure 4, is 42.9 inches (109 cm) long and 19.1 inches 
(48.5 cm) wide. It has seven alternating up and down semicircular 
beads separated by about 0.4 inches (1 cm) wide flats. The perimeter 
of the panel is flat to permit mounting to the spar daps and rib caps 
of a wing. 

Doublers made from Rend 41 sheet stock were spot welded to the 
ends of the panel on both sides, tripling the nominal thickness of the 
ends. These doublers extend about 10.7 inches (27.2 cm) towards the 
center of the panel, gradually reducing in thickness as they progress 
down the flats. The doublers reduce the possibility of local 'end 
failures and help to distribute the load more uniformly into the panel. 

Provisions were made to attach heat shields to the panels at 
eight locations two of which are pointed out in figure 4. The remain- 
ing six attachment points are symmetric to those shown. 

The beaded panel dimensions were derived using a computerized 
optimization program (ref. 2). The optimization program varied such 
parameters as panel length, width, number of beads and thickness to 
derive a least weight configuration capable of carrying prescribed 
mechanical and thermal loads (based upon semiclassical analysis). 
Strength interaction curves such as that shown in figure 5, were then 
made for each panel configuration. The various combinations of 
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Figure 4. The hypersonic beaded skin panel used for this project. 
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Rend 41 Material and Formabilicy Properties 


Tensile specimens, shown in figure 6, were supplied with the beaded 
panels. These specimens and the panels were cut from the same sheet 
stock. A number of them were cut in the direction of rolling; an 
equal number were cut perpendicular to that direction. Six specimens, 
three of each type, were tested in a universal testing machine to ex- 
perimentally determine the modulus of elasticity and the 0.2% offset 
yield strength of the Rend 41. Figure 7 is a typical stress-strain plot 
using test data. The average modulus of elasticity was found to be 
30.422xl0 6 psi (2.10xl0 11 N/m 2 ) compared to 31 . 6x1 0 6 psi (2.18xlO^N/m 2 ) 



Figure 6 . Rene 4l tensile specimen dimensions 0.037 
inches (0.094 cm) thick. Dimensions in 
inches and (cm). 


as reported in reference 1 A . Table 1 summarizes the results of the 
six tensile specimen tests In table 2, a number of material properties 
are summarized as reported in reference 14. 

The beaded panels were formed in a 5 million pound (2.22x10 N) 
hydraulic press. At least two and sometimes four anneals were required 
before fully developing the bead (ref. 15). The stretch forming process 
reduced the thickness of the bead from an original thickness of 0.037 
inches (0.094cm) to about 0.028 inches (0.071 cm). 

Due to the extreme hardness of Rend 41, standard high speed steel 
drills could not be used. Strict drilling procedures in addition to 
cobalt drills had to be used to prevent work hardening the Rend 41 and 
to obtain maximum life from the drill bits (ref. 15). 




0 200 400 600 800 1000 X 10 1 

Strain, pm/in or pern/ cm 


Figure 7. Typical stress-strain curve for 
Bene 41 tensile specimen test. 
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Table 1. Experimentally determined modulus and 0.2 % offset yield 


strain of Rend 41. 


Specimen 

ID 

Modulus of 
Elasticity 

0.2% offset 
Yield strain 

lbf/in^ (N/cm^) 

vnn/in (ycm/cm) 

003-L-8* ' 

30.205X10 6 
{20.826x1 0 6 ) 

6,350 

003-L-9 

30. 358x1 0 6 
(20.931x1 0 6 ) 

6,450 

003-L-10 

30. 091x1 0 6 
(20.747x1 0 6 ) 

6,200 

Q03-T-8** 

30. 668x1 0 6 
(21. 145x1 0 6 ) 

6,250 

003-T-9 

30. 126x1 O 6 
(20.771xl0 6 ) 

6,600 

003-T-10 

31. 082x1 0 6 
(2 1. 430x1 0 6 ) 

6,400 

Average 

30. *22x1 0 6 
(20.975X10 6 ) 

6,380 


*L-Longitudmal specimens **T-Transverse specimens 
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Table 2. Material properties of Rend 41 as reported m reference 14 


Property 

Val ue 

Modulus v of elasticity (Tension), E 

31. 6x1 0 6 lbf/in 2 
(21. 8x1 0 6 N/cm 2 ) 

Modulus of elasticity (Compression), E c 

31. 6x1 0 6 lbf/in 2 
(21. 3x1 0 6 N/cm 2 ) 

Shear modulus, G 

12.1xl0 6 lbf/in 2 

» 

(8.3xl0 6 N/cm 2 ) 

Density, p 

.298 lbm/in 3 
(.008 kg/cm 3 ) 

Tensile ultimate strength 

185,000 lbf/in 2 
(127,550 N/cm 2 ) 

Tensile yield strength 

132,000 lbf/in 2 
(91,000 N/cm 2 ) 

Compressive yield strength 

< 

141,000 lbf/fn 2 
(97,200 N/cm 2 ) 

Poisson's ratio, v 

.31 






















PRETEST PREPARATION 


Beaded Panel Measurements 

Detailed drawings accompanying the beaded panels to NASA specified 
the various dimensions of the panel which are summarized in xable 3. 
Also included in table 3 are the values obtained from direct measure- 
ment which vary somewhat from the specifications. The measured values 
were used in all analyses of this project. 

The buckling characteristics of any structure are affected by 
eccentricities. Therefore, measurements were made at over 100 loca- 
tions on the beaded panel to determine the magnitude of the manufactur- 
ing eccentricities present in the panel. The edge stiffeners (which 
will be discussed in a later section) were attached to the panel during 
these measurements. The maximum out-of-olane eccentricity of the panel 
was found to be only 0.019 inches (0.048 m). 
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Edge Stiffeners 


An attempt was made in this project to stiffen the edges of the 
panel in a manner which would closely approximate the stiffness con- 
ditions of a wing mounted condition. Therefore, the edge stiffeners 
which run parallel to the bead as shown in figure 8 were sized with 
the intention of not only preventing local edge failures but also 
simulating the stiffness of adjacent spars, and panels in the HWTS. 
These stiffeners were made in the shape of Z-sections from annealed 
stainless steel and mounted on the heat shield side of the panel as 
shown in figure 9. 

The Z-sections were designed with the aid of Nastran, a finite 
element structural analysis program. Nastran was used in an iterative 
manner utilizing two finite element models entitled EDGE! and EDGE2. 

Model EDGE!, shown in figure 10, consisted of a quarter Dane!, 


K 

2.04 

(5.18) 


r~ 

< 42 


(108.5) 






1.42 

(3.61) 


*1 


K-0.85 

(2.16) 


Figure 8. Dimensions of the Z- section edge stiffeners made from 

annealed stainless steel. Dimensions vn inches and (cm). 
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• Side stiffeners 3 end supports and stabilizing rod mounted to the beaded panel . 




Figure 9 




spar cap, spar web and an adjacent panel assembly. This model was 
intended to be a representative section of an actual hypersonic wing. 

Model EDGE2, as shown in figure 11, consisted of the same quarter 
panel as used in EDGE! but bar elements replaced the spar cap, spar 
web and adjacent half panel assembly. The design procedure employed 
was to apply identical compressive loads (parallel to the beads) to 
both models, then adjust the sizes of the bar elements until the out- 
of-plane displacements of model EDGE2 were comparable to model EDGE!. 
This procedure provided the dimensions of bar elements which approached 
the bending stiffness of the spar cap, spar web and adjacent half 
panel assembly. 

The results of this iterative procedure are shown in figure 12. 
The curves, s,hown in figure 12, represent about twenty iterations 
and are the best possible correlations obtainable. The Z-sections 
were dimensioned on the basis of the computer run of model EDGE2 cor- 
responding to the curve shown in figure 12. 




Figure 1 1 . 


Nastvan model EDGE2. 
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End Supports 


The end supports shown mounted on the panel (previously shown 
in figure 9) were made from machined flat tool steel bars. The bars 
had a rectanqular cross section of 1 inch by 0.5 inch (2.54x1.27 cm) 
and were flat to within 0.001 inches/inch. When mounted on the panel 
the end stiffeners served two functions. They provided a surface 
approximately 1.10 inches (2.79 cm) wide (including stiffe'ners and oanel 
thickness) through which the load was transferred into the panel and 
eliminated warping of the end of the panel. 

After the side and end stiffeners were mounted on the panel, the 
entire assembly was mounted in a milling machine square with the 
cutting tool. The end stiffeners were then milled off parallel with 
one another and oerpendicular to the beads. This orocess was necessary 
to ensure that bending loads would not be introduced into the oanel 
due to misaligned ends. The side and end stiffeners were not removed 
after this process had been completed. 



Strain Gage Instrumentation 


Figures 13 (a) and (b) show the location of the 39 strain gages 
mounted on either side of the panel to measure strains. The gages are 
mounted on cross sectional lines corresponding to 1/4, 1/3, and 1/2 
the panel length. Of the 39 gages, 33 were standard axial gages. The 
remaining six gages were grouped into threes and used as equiangular 
rosettes [labeled 34 and 37 in figure 13(b)]. All of the gages were 
attached using standard strain gage adhesives. 

Positioning of Strain Gages on Beads 

As discussed in the Semiclassical Buckling Analysis section (later 
in this paper), a diagonal mode of local instability which occurs between 
two adjacent beads has been suggested as a possible failure mode (ref. 2, 
3, 4). For the beaded panels, the buckling load for this mode was de- 
termined for a value of 0 2 equal to 12.875° (see fig. 14). Thus the 
majority of gages mounted on the beads were mounted at about 12. 8 9 off 
of the bead peaks as shown in figure 14. This placed the gages at 
locations that would optimize their sensitivity to both the proposed 
diagonal and general instability modes. 

The strain gages were mounted at three cross sections of the panel 
corresponding to 1/2, 1/3, and 1/4 panel lengths. The majority of the 
gages were mounted on the 1/2 panel cross section, in anticipation of 
maximum panel deflection at that location. At the 1/3 and 1/4 panel 
cross sections the gages were clustered around the center three beads 
[see fig. 13(a) and (b)]. 
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TaJ Side 7. 



(ZU Side 2. 


Figure 73. Strain gage instrumentation locations. Note: Strain 

gage identification numbers should be preceeded by the 
nimber 4 to correspond to the remainder of the report. 
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Figure 14. Location of strain gages mounted on the beads for maximum 
sensivity to the diagonal mode of instability proposed in 
the semialassioal analysis. Note: 02 = 12.8°. 
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Displacement Transducer Instrumentation 


Displacement transducers (DT's) were located on panel cross 
sections at 1/4, 1/3 and 1/2 panel lengths as were the strain gages. 

As shown in figure 15, five DT's were attached at the 1/2 panel cross 
section in anticipation of maximum out of plane displacement at that 
location. Three additional DT's were placed at the 1/4 panel cross 
section and one DT was located at the 1/3 panel cross section. All 
nine of these DT's were positioned to measure displacements perpendic- 
ular to the plane of the beaded panel. A tenth DT was used to measure 
longitudinal compression of the panel. 

Moir6 Fringe Technique Preparations 

The photographic material to which the Moir£ fringe grid lines 
were applied was only 0.007 inches (.018 cm) thick (see appendix C). 

The plastic was attached to a 0.25x20x42 inch (.62x50.8x106.7 cm) 
sheet of plate glass to enable mounting the grid plane in front of the 
panel. Mineral oil was used as an adhesive between the glass and the 
plastic. Excess oil was squeezed from between the glass and the plas- 
tic creating a thin uniform adhesive bond. 

The Moir6 fringe glass was supported by aluminum bars which ran 
the length of the glass. A 0.25 inch (0.64 cm) groove was cut into 
each of the aluminum bars and the glass fitted and glued with silicone 
rubber cement into the grooves. The aluminum bars were then attached 
to the panel by aluminum brackets such that the grid was maintained 
parallel and at a fixed distance from the ends of the panel. The bottom 
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Figure 15. Displacement transducer (DT) locations. 
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brackets were rigidly attached to the panel and the top ones were 
free to slide, thus preventing the glass from taking any load and 
possibly breaking, or from deflecting and causing errors in the measure- 
ments. The distance from the top of a bead to the surface of the glass 
was approximately 0.25 inches (0.64 cm) in the unloaded condition. 

Prior to mounting the glass to the panel, the panel was painted white 
to create a greater contrast between the dark and light fringes. 

The camera film plane was located about 60 inches (152.4 cm) 
from the surface of the glass (S^ in fig. Cl in appendix C). The dis- 
tance between the camera and the light source, was 60 inches (152.4 cm) 
(S2 in fig. C2 in appendix C). 

According to reference 12, it is necessary to use a point source 
of light when the field of interest is large. Therefore, a photo- 
graphic flash with a 1 inch by 0.2 inch (2.54x. 51 cm) iris was used 
as a light source for this project. The iris effectively created the 
necessary point source of light. Figure 16 is a photograph of the 
entire photographic system in place. 

Test Equipment 

The panel was tested in a universal compression-tension testing 
machine. Figure 17 shows the panel mounted in the machine. The 
platens which come into direct contact with the panel were specially 
made and machined flat to within 0.001 inches/in. The platens helped 
to insure that the load introduced into the panel was purely axial in 
nature and that bending loads due to misaligned heads would be elimi- 
nated. The bottom platen rested on a spherical seat which insured 




Figure 17. Beaded panel test assembly ^ 

installed in the testing machine. 




Test Procedure 


A total of 19 tests were oerformed on the panel up to the failure 
test of 48,600 pounds {216,184 N). The first two tests were system 
checkout tests of 2,000 and 10,000 pounds (8,900 and 44,480 N) respec- 
tively. The remaining 17 tests were buildup tests to failure. 

A typical test would begin by warming up the testing machine 
for 30 minutes. Before loads were applied to the panel, Moird fringe 
calibration photographs ( were taken. In addition, strain gage and dis- 
placement transduder zeros were recorded by the data acquisition 
system. After data sampling was started at prescribed rates, loads 
were applied to the panel. At predetermined load points, the load 
was held constant and Moird fringe photographs were taken. This pro- 
cess was repeated until the maximum load was reached. 
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DISCUSSION OF EXPERIMENTAL RESULTS 


Moird Fringe and Displacement Transducer Results 

The Moire fringe photographs proved to be extremely useful for 
identifying the elastic buckling load and the mode shaDe. The Moird 
fringe photographs show subtle changes in panel curvature which might 
not have been otherwise detected. Experimental data, including Moird 
fringe data, is included in appendix A. 

Figure Al(a) is a calibration photo taken while no load was being 
applied to the panel. Calibration bars (with different slopes) are shown 
in the upper and lower right corners of the nanel between the glass and 
the panel. The bars were used to verify uniform calibration from end to end. 
Note that the beads and flats are all straight in -figure Al(a) as indicated 
by constant fringe patterns on bead peaks and flats. 

Figure Al(b) shows the panel under an applied load of 24,000 pounds 
(106,760 N). A very slight curvature of the center two flats is visi- 
ble, compressing the sides of the center bead. This inward deflection 
is even more apparent at 36,000 pounds (160,140N), as shown in figure 

i 

Al(c). The sides of the other beads are similarly compressed inward or 
spread outward but to a lessor extent than the center bead. This deflec- 
tion represents lateral distortion of the panel across section due to out- 
of plane bending. 

The first visible indication of elastic buckling of the panel 
occurred at 4 q, 000 pounds (177,930 N), as shown in figure Al(d). The 
flat immediately to the left of the center bead has changed its direction 
of lateral deflection and is now moving outward, away from the center 
bead, instead of its original inward direction. At 42,000 pounds 
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(186,825 N) , shown in figure Al(e), this change in curvature is distinct, 
while at 44,000 pounds (1 95,520 N), shown in figure Al(f), it is quite 
pronounced This mode of buckling is similar to the diagonal mode of 
instability discussed later in the semiclassical analysis section 

In figure Al(g) the panel is under an applied load of 48,500 pounds 
(215,740 N) . Severe curvature of the center bead and the inner most flats 
is clearly visible. Unlike figure Al(f) which shows the panel under an 
applied load of 44,000 pounds (195,720 N) , the beads immediately adjacent 
to the center bead are beginning to exhibit curvature similar to that 
of the center bead. The remaining beads, however, remain relatively 
straight. This curvature of the panel center portion suggests that the 
center portion of the panel is carrying less of the applied load, having 
transferred some of the load to the outer portions of the panel. This 
load transfer was further substantiated by the strain gage results 
which will be discussed in the next section. 

Ultimate panel failure occurred at a load of 48,600 Dounds (216,184N). 
Figure Al(h) is a Moire fringe photo taken after failure. Figures Al(i) 
and (j) are photos of the panel after failure. Permanent deformation 
is visible in these photos. These figures show that the panel suffered 
catastrophic local failure at the center of the panel. Local failures 
are also visible between the fasteners on the panel edges. 

Out-of-plane displacement measurements using the Moird fringe tech- 
nique were made at 49 locations as shown in figure A2(a). All of the 
measurement points were located on the peaks of the seven beads. Note 
that the panel deflected away from the Moi rd fringe glass. In figure 
A2(b) through (h), plots of out-of -plane displacements are shown for 
each individual bead. Where possible, displacement transducer data is 
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also plotted. The curvature of these plots indicate that the panel was 
loaded eccentrically which caused out-of-plane bending of the panel 
even at small loads. The eccentric behavior was the result of loading 
the panel through the flat end of the panel , rather than the neutral 
axis of the panel/side stiffener assembly. The maximum out-of-plane 
displacement at a load of 48,500 pounds (215,740 N) was 0.432 inches 
(1.097 cm) as shown in figure A2(e). This displacement occurred on the 
center bead at point 25 in figure A2(a). 

Strain Rage Results 

Initial interest in the strain gage results was focused upon the 
degree to which uniform loading had been accomplished. Uniform loading, 
in this case, refers to a uniform load across the entire width of the 
panel. The two rosettes (gages 434 and 437) mounted on the panel flats 
provided part of this information by making it possible to resolve the 
axial strains into principle strains and principle directions. At all 
load levels the difference between the measured longitudinal strain and 
the calculated principle strain was negligible. In addition, the direc- 
tion of the principle strains varied by a maximum of only 5° from an 
axis parallel to a bead. This small variation indicated that the load 
was introduced into the panel in a uniform manner and that there was 
virtually no shear. Sample data is shown in table 4 for a load of 
20,000 pounds (88,960 N). 

Another indication of the uniformity of the load is illustrated 
by the data shown in figure 18. This figure shows the strains from 
all strain gages recorded at a load of 2,000 pounds (8,900 N). At this 
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low load very little bending is present, thus the indicated strain is 
due primarily to axial compression. Similar comparisons of strains on 
beads at higher loads cannot be made due to the increased effects of 
bending. However, strain comparisons can be made of the responses of 
gages mounted on flats at higher loads. Figure 19 shows the strain 
measurements made at four load levels on the flats at the center cross 
section. The maximum difference between any two gages at a particular 
load is 150 microinches/inch. The data in this figure shows that (1) 
the strains are increasing in nearly equal increments with each load 
level and (2) that the compressive load is uniform across the beaded 
panel cross section up to the onset of elastic buckling (about 40,000 
pounds or 177,930 N). 


Table Principle strains and their directions at a load of 20,000 


pounds (88,960 N). Strains in pinches/inch. 


Rosette 34 

Rosette 37 

measured 

calculated 

measured 

calculated 

Leg A -876 

e-j -876 

Leg A -864 

e-j -868 

Leg B - 25 

e 2 " 252 

Leq B - 78 

e 2 -272 

Leg C - 35 

Y - 12 

Leg C + 48 

Y H5 

princiole angle 89.7° 

principle angle 86° 
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Figure 19. Strain measurements recorded by gages mounted on the 
center section of the beaded -panel. 
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Individual strain gage plots are shown in figures A3(a) through 
(w). Most of these plots are nonlinear. These non! mean ties are due 
to out-of-plane bending which is the direct result of eccentric loading 
of the panel . 

This nonlinear load deflection response is typical for most column 
structures under compressive load (ref. 17) and is the reason for 
the difficulty associated with Pinpointing elastic buckling loads from 
strain gage plots alone. A column with no eccentricities and concentric 
loading would have a load-bending deflection response similar to that 
shown in figure 20(a). This figure represents a perfect column under 
compressive load, where the column simply compresses until the buckling 
load is reached. For a column with eccentricities (or eccentric loading) 
the load-bending deflection response is represented by figure 20(b) 
where a. is a measure of the eccentricity and a,<ct <a etc. Since the 

i 12 3 

beaded panel was eccentrically loaded, the strain gage plots are similar 
to figure 20(b). In those cases where the strain gage measured axial 
compression plus compressive bending, the resulting plot is similar to 
that shown in figure A3,(b). For the case where the. gage measured axial 




Deflection Deflection 

Figure 20 . Load bendvng-de flection responses for structures with (a) 

no eccentricities and (b) with eccentricities of magnitude cu. 



compression plus tensile bending, the plots are similar to that shown 
in figure A3(h), (note the reverse curvature). A precise determination 
of the buckling load for such plots is difficult because the curves do 
not exhibit a pronounced change which identifies buckling. In the case 
of the beaded panels, the situation is further complicated by the fact 
that the elastic buckling mode, as shown in the Hoird fringe photos, has 
a lateral component which is perpendicular to the deflection due to the 
eccentric loading. The best indication of buckling from the strain gage plots 
is given by gages which are mounted on the flats which are relatively 
insensitive to out-of-plane deflection (bending), since they are much 
closer to the neutral axis of the panel (gages 426 and 428 in figures 
A3(q) and (r), for example). 

Figure 21 is a plot of the average strains recorded by all of the 
gages mounted on the flats (between beads) at the center cross section 
of the panel (gages 425, 426, 428, 429, 434, and 437). Up to about 
40,000 pounds (177,930 N) the average strain gage response is linear. 

This load corresponds to that at which elastic buckling of the center 
of the panel occurred as shown in the Moiri fringe photos. The average 
strain at 40,000 pounds (177,930 N) from figure 21 is 1600 microinches/ 
inch and does not increase appreciatively at loads above 40,000 pounds 
(177,930 N). 
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Force/Stiffness Results 


Force/stiffness plots are shown in figures A4(a) through (j). 

With the exception of plots (c) and (j), the curves were extrapolated 
to indicated failure loads using a second order Lagrange polynomial 
(ref. 18). The accuracy of the predictions are affected by the extra- 

i 

i 

polation procedure used i.e. , linear, second order, third order, etc. 

A second order procedure was used for this project since the curves 
extrapolated were generally quadratic in nature. 

The most important results of this analysis are shown in figures 
A4(c) and (j). In these two plots, very pronounced inflections occur 
at about 40,000 pounds (177,930 N) as indicated in the figures. 
According to reference 11 these inflections indicate changes in the 
mode of deflection. More specifically, the inflections are caused by 

4 

elastic buckling of the panel which occurred when the flats on either 
side of the center bead began moving in the same lateral direction. 
Since the inflection points themselves indicate elastic buckling, the 
curves were not extrapolated. It should be pointed out that the 
significance of these inflections in figures A4(c) and (j) was re- 
alized only after correlations between Moird fringe and strain gage 
data were made. 

Other results of particular interest are the plots of gages at the 
center cross section of the panel. Extrapolations of plots A4(a), (b), 
(d), (e), and (f) all intersected the load axis at between 50,850 to 
55,000 pounds (226,200 to 240,650 N). These indicated buckling loads 
are for the mode of deflection associated with out-of-Dlane deformation 



only. Had the panel not failed elastically at 40 s 000 pounds (177,930 N) 
in a mode characterized not only by out-of-plane deflection, but lateral 
deflection also, it may have supported loads up to the indicated buckling 
load. Force/stiffness predictions from gages mounted at the 1/4 and l/3„ 
panel lengths ranged from 53,750 to 61,750 pounds (239,100 to 274,680 N). 
The larger predictions came from gage pairs mounted at locations which 
were relatively unaffected by the elastic failure which occurred at the 
center of the panel. These results point out that the strain gage pairs 
must be in close proximity of the failure for the most accurate results. 

Two points should be made regarding similarities between the experi- 
mental results, i.e., modes of deflection, and the modes of instability 
suggested by the semi classical analysis (discussed later in this paper). 
One suggested mode of instability used to analyze the panels, is charac- 
terized by a lateral deflection corresponding to the lateral deflection 
which occurred beginning at 40,000 pounds (177,930 N). In the semiclas- 
sical analysis this mode of instability has been called the diagonal mode. 
Secondly, the out-of-plane deformation of the panel corresponds to one of 
the suggested instability modes in the semiclassical analysis, known as 

the general instability mode. The semiclassical analysis section con- 

> 

tains a further discussion of the ( suggested modes of instability. 



Edge Stiffener Performance 


The edge stiffeners (Z-sections) used to support the edges of the 
beaded panel were intended to simulate the stiffness of a wing mounted 
condition. The out-of-plane displacements of the center of the beaded 
panel edge were compared with similar measurements taken from a panel 
mounted in the hypersonic wing test structure. The deflection of the 
beaded panel /Z-section edges were larger than the corresponding HWTS 
panel edges. Thus, the Z-sections were stiffened with 1x1x42 inch 
(2.54x2.54x106.68 cm) steel bars which were bolted to the free edqes 
of the Z-sections as shown in figure 22. The out-of-plane displacement 
of the beaded panel /stiffened Z-section edge (at the center) for a com- 
pressive load of 950 Ibf/in (1660 N/cm) was 0.015 inches (0.038 cm). 

The corresponding measurement of the beaded panel mounted in the HWTS 
was 0.024 inches (0.061 cm). 

The reason for the conservative panel /Z-section deflection are two- 
fold. First, an incorrect load level was taken from the HWTS data and 
used for comparison purposes. Secondly, the Nastran calculations of 
out-of-plane deflections were found to be inadequate (see section on 
computer analysis). However, the edge stiffeners did prevent premature 
edge failures and in general performed satisfactorily. 

It was discovered during the initial tests of the beaded panel, 
that the free edges of the Z-sections tended to pull away from one 
another (buckle laterally) under load. To prevent a premature lateral 
buckling failure of the Z-sections, 0.5x0.25 inch (1.27x0.64 cm) steel 
straps were used to tie the Z-secticns together, as shown in figure 22. 


45 



original mo® ® 
OF POOR QUALrTY 



Figure 22. Additional edge support bars and straps 
mounted on the panel. 
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NASTRAN COMPUTER MODEL DESCRIPTIONS AND DISCUSSION OF RESULTS 

Introduction 

Nastran structural analysis (ref. 13) was used as a tool to com- 
pute deflections, stresses and the buckling strength of Dortions of 
the panel as well as the entire panel. Secondly, Nastran was used as 
an aid to design the edge stiffeners as previously discussed. 

A total of five structural models were made and used in the analy- 
sis process. Models EDGE2 and EDGE3 were 1/4 and 1/2 panel models, 
respectively. Each of these models were used to determine the strength, 
deflection and buckling characteristics of the beaded panels. The 1/2 
panel model EDGE3 was made after the results of the buckling analysis 
done with the 1/4 panel model EDGE2 were found to be inadequate. The 
remaining three models entitled BEAD, FLAT and DIAG, were used to ana- 
lyze instability modes suggested in the semiclassical analysis classi- 
fied as bead, flat and diagonal modes of instability, respectively. 

The results of the computer analysis done with models EDGE2 and 
EDGE3 will be compared with the experimental results previously dis- 
cussed. The other three models, BEAD, FLAT and DIAG will be compared 
to the semiclassical analysis only. 
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Finite Element Model EDGE2i Description 
and Results of Analysis 

Figure 23 is a computer generated clot of model EDGE2. It is a 
quarter panel model, and takes advantage of the two lines of symmetry 
of the beaded panel. It consisted of 306 elements, 298 grid points and 
1443 degrees of freedom. EDGE2 utilized dimensions and thicknesses of 
the panel as determined by direct measurements. In addition, the eccen- 
tricities which were measured after the edge stiffeners and end supports 
had been mounted were also incorporated into the model. 

' Figure 23 also shows the direction of the loads applied to model 
EDGE2. The letters along the sides of the model indicate the restricted 
degrees of freedom. For example; x indicates that translation in the 
X-direction is restricted, Rx indicates rotation about the X-axis is 
restricted. The boundary conditions along the right and lower sides of 
the panel were relatively straight forward to define, since the panel 
attaches to spar and rib caps at these boundaries. However, the boundary 
conditions along the cut edqes (lines of symmetry) were not as easy to 
formulate. This was especially true of the instability or buckling 
analysis. The mode shape prior to testing was expected to be the gen- 
eral instability mode (See semi classical analysis section). Therefore, 
the boundary conditions along the cut edges of the panel model were 
varied in a trial and error Drocedure until olots of the buckled shape 
resembled the general instability mode. ' * 

The results of the buckling, analysis performed by model EDGE2 were 
poor. An elastic buckling load of 31,700 pounds (141,000 M) was calcu- 
lated which compares with an actual elastic buckling load of about 


49 




Figure 23. Finite element buckling model EDGE 2. The restricted degrees 
of freedom indicated along the sides of the model are for 
buckling analysis only. 
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40,000 lbs (177,930 N). The predicted load was 26% lower than the 
actual failure load. This large descrepancy prompted the development 
of the 1/2 panel model, EDGE3. The reason for the large difference be- 
tween the results of the buckling analysis utilizing model EDGE2 and 
the experimental results is very probably caused by the problems associ 
ated with prescribing boundary conditions along the cut edge of the 
model (i.e., the line of symmetry parallel to the center bead). 

It should be noted that initially the analysis done with model 
EDGE2 was done assuming a perfectly flat structure, i.e., with no eccen 
tricities. However, the difference in results between the analysis 
including eccentricities and the analysis without eccentricities was 
insignificant. 

Finite Element Model EDGE3* Description 
and Results of Analysis 

Model EDGE3 shown in figure 24 was a full half panel model con- 
sisting of 920 elements, 842 grid points and 4591 degrees of freedom. 

As m model EDGE2, the dimensions and thicknesses used in model EDGE3 
were measured directly from the panel. Measured eccentricities were 
not included since previous experience with model EDGE2 had shown that 
inclusion of the eccentricities had an insignificant effect on the 
results of the analysis. 

Figure 24 also shows the direction of the loads which were applied 
to the panel and the restricted degrees of freedom. The same degrees 
of freedom used in model EDGE2 were also used in EDGE3. 

The buckling analysis using EDGE3 gave an elastic buckling load 




Figure 24. Finite element buckling model EDGE3. The restricted degrees of freedom indicated along 
the sides of the model are for buckling analysis only. 
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of 50,000 pounds (222,400 N). This calculated load compared to an ac- 
tual elastic buckling load of about 40,000 pounds (177,930 N) repre- 
senting an error of 25%. The magnitude of the error reflects the accu- 
racy of the model for the buckling analysis. 

The Nastran static analysi s program exhibited an inability to account 
for the nonlinear effects of out-of-plane bending with loads applied in 
the plane of the panel. Thus comparisons of calculated and experimental 
stresses across the entire panel, could only be made at relatively low 
loads, where the effects of bending were small. In addition, stress com- 
parisons at higher loads could only be made on the flats where the effects 
of bending were not predominant. For example, figure 25(a) is a plot of 
the stresses at the center cross section of the panel, at a load of 
10,000 pounds (44,480 N) using models EDGE2 and EDGE3 as well as the 
experimentally derived stresses. At 10,000 pounds (44,480 N) a good cor- 
respondence exists between the analysis and experimental results across 
the entire panel . However, the experimentally derived stresses on the beads 
of the panel are noticeably affected by out-of-plane bending whereas the 
computer analysis results are not. Figures 25(b) and (c) are plots of the 
stresses at the center cross section at 24,000 and 36,000 pounds (106,760 
and 160,140 N), respectively. Pronounced bending effects can be seen in 
the experimental results; these effects are not accounted for by Nastran 
static analysis (with applied in-plane loads). On the flats, where bending 
has less effect, good comparisons are possible even at the higher loads. 

The results of the analysis using these finite element models point 
out the need for further research in the area of finite element buckling 
analysis for the beaded panels. The development of a full panel model 
would be one possibility for further study. 
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Figure 25. Comparisons of static stresses computed by finite element models EDGE2 and EDGE3 with experimental 

results. 
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Finite Element Models BEAD* FLAT and DIAG: 

Descriptions and Results of Analysis. 

The remaining three Nastran structural models, BEAD, FLAT and 
DIAG, were models of an individual panel bead, a flat between beads, 
and of a cross section from approximately peak to peak of two adjacent 
beads, respectively. These models were used only for comparison pur- 
poses with semi classical analysis. 

The finite element model BEAD is shown in figure 26. It con- 
sisted of 447 elements, 497 grid points and 1347 degrees of freedom. 
Boundary conditions for BEAD were prescribed to simulate simply sup- 
ported edges. The restricted degrees of freedom and direction of the 
applied loads are shown in figure 26. 

The Nastran model FLAT was a model of the flat between two beads. 
It was a very simple model consisting of 84 elements, 170 grid points 
and 344 degrees of freedom. The boundary conditions used were identical 
to those shown along the edges of the model BEAD shown in figure 26 and 
represented simply supported edges. 

The last of the three section models was DIAG which is shown in 
figure 27. It consisted of 320 elements, 374 grid points and 1694 
degrees of freedom. As in the case of models BEAD and FLAT, the bound- 
ary conditions on the edges of DIAG represented simply supported edges. 
Each of the foregoing models was loaded in compression parallel to thei 
long axis. 

The results of the analysis using these three section models were 
compared to semiclassical results only, (which will be discussed later) 
since no individual section tests were performed. The comparisons are 
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Figure 26. Finite element buckling model BEAD, 
simply supported edges. 
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Figure 27. Finite element buckling model FI AG. The restricted degrees of freedom indicated represent 
simply supported edges . 
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summarized in table 5 The best comparison occurs with the flat in- 
stability modewith a difference of 16%. The bead and diagonal modes 
exhibit differences of 23% and 58%, respectively. The reason for these 
poor comparisons is not known. 

Table 5. Comparison of section model buckling results with semi- 
classical analysis. 


Instabil ity 
mode 

Seimclassical 
buckling load 

Nastran model 
buckling load 

Percent 

difference 

Bead 

343,200 lbf/in 2 

264,600 lbf/in 2 ' 

23 


(236,600 N/cn 2 ) 

(182,400 N/cm 2 ) 


Flat 

687,000 lbf/in 2 

579,700 lbf/in 2 

16 


(473,600 M/cm 2 ) 

(3 Q 9,7QQ N/cm 2 ) 


Diagonal 

42,700 Ibf 

67,400 lbf 

58 


(190,000 N) 

(299,700 N) 
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SEMI CLASS I CAL ANALYSIS 


Introduction 

The purpose of this section is to present a summary of the semi- 

classical analysis of the beaded panel buckling characteristics given 

* 

in references 2, 3, and 4. The four beaded panel buckling modes which 
were identified in the references include general instability, flat in- 
stability, bead instability and diagonal instability. The diagonal 

mode of instability is a localized instability between the peaks of 

$ 

adjacent beads. In all four instances, semiclassical buckling theory 
assuming simply supported edges was used in the analysis. The use of 
this simplified approach on a problem with the complexity of a beaded 
panel is questionable. 

Since this project is concerned solely with compressive loads on 
the panels just the pertinent equations will be presented. A more 
detailed analysis may be found in references 2, 3 and 4. 

Section Properties of the Panel 

Prior to examining the four instability modes of the panel, it is 
necessary to define a number of section properties to be used in the 
various buckling equations. The first of these properties are the 
flexural rigidities of the beaded panel with respect to the X, Y and 
£ axis shown in figure 28. The equations of flexural rigidity for 
bending moments along the X and Y axes and twisting of the XY plane 
are: 
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1.045 in 
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a 


X 


c 


= .8190 
(2.080 


in 
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t„ = .0345 in 
(.0876 cm) 


= 1.239 in 
(3.147 cm) 

t = .0285 in 
(.0724 cm) 

01 = 77.5° 



Figure 28. Section ’parameters of the beaded panel 
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The diagonal mode of instability occurs at a critical value of 
02 » shown in figure 28, where the area moment of inertia of the diagonal 
section is minimum. The flexural rigidity equations about the £Y plane 
for buckling across the diagonal defined by angle 0 2 are: 

D I ~ 



n = jL. 
11 S 6 

Et 3 

(3) 

12(1 - v 2 ) 

„ ,fi 

Et 3 



12 " 5 12(1 + v) 


in equation 3 includes only the material within the dimensions 
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of 6 and is defined by: 


= 1 (I xx Cos 2 0 3 +.I ZZ Sin 2 0 3 - I Y7 Sin 26 3 ) 


l XZ 


(4) 


A /V 


and I xx , I zz , and I xz are defined as: 


I y Y = al v - tFT{(0.5 + Cos^ei)e 2 - Sin 6 2 (2 Cos Qi-0.5 Cos 9 2 )} 

AA XX 


l TL = tR *(f + a2 )( e i- 0 2) ~ V (Sm 2G!-Sin 20 2 ) + 


+ 2aR{Cos 6i - Cos 0 2 ) 


t F b 3 
24 


(5) 


i xz = tR2{a ( Sln 0 i- sin 0 *) - f ( s1n2 01 " Sir * 2 02 ^ 


- R Cos 0i (Cos 0i - Cos 0 2 ) - (0i - 0 2 )a Cos 6i> 
Numerical values for equations 1 through 5 are summarized in table 6. 
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Table 6. Numerical values of the 


Section 

Property 

Numerical 

Value 

D i 

49.286 lbf-in 
(5.569 N-m) 

°2 

33.75x1 0 4 lbf-in 
(38.1 3x1 0 4 MN*m) 

°3 

59102 lbf -in 
(6.67 N-m) 

X 
X 
» — 1 

1.1 09x1 0~ 2 in 4 /m 
(.462 cm 4 / cm) 

1 zz 

1.1 64x1 0‘ 2 in 4 / in 
(.484 cm 4 / cm) 

/s. 

'xz 

1. 029x1 0“ 2 in 4 /in 
(.428 cm 4 / cm) 


id panel cross section parameters 


Section 

Property 

Numerical 
Val ue 

D i 

5,415 lbf-in 
(611.8 N-m) 

D n 

59.51 lbf-in 
(6.72 N-m) 

°12 

48.88 lbf-in 
(5.52 N-m) 

i /*> 

I 

XX 

9.256x1 0‘ 3 in 4 / in 
(.039 cm 4 /cm) 

A 

1. 780x1 0" 4 i n 4 /i n 
(7. 409x1 0~ 3 cm 4 / cm) 





























Semiclassical Buckling Analysis for Pure Compression 


Stability equations for the beaded panel may now be summarized. 

It is important to reiterate that the equations which follow are based 
upon simplified classical theory. 

General instability is analyzed by assuming that the entire panel 
* is a simply supported wide column under compressive load. From refer- 
ence 2 the critical axial buckling load is given by: 


N 


YCr 



' ( 6 ) 


Diagonal instability, a local instability (of the simply supported 
panel cross section) between the oeaks of adjacent beads is given as 
follows: 


YQCr 


2^2 ^ D I D II^ + D 12* s^ s 


(7) 


The critical diaqonal buckling load occurs when the diagonal 
cross section has a minimum moment of inertia. The angle 02. in figure 
28 which defines the boundaries of the diagonal cross section, was 
varied from 0 to 77.5 degrees m small increments and the critical 
buckling load calculated. The results are shown in figure 29 which 
is a plot of the diagonal buckling load versus angle 0 2 . Given the 
geometry of the beaded panel used in this oroj'ect, the angle 0 2 is 
about 12.875 degrees from the peak. 

Instability of the flats between the beads is the third mode to 
be considered. From reference 2, the compressive buckling Formula 



for flat instability, assuming simply supported edges, is given by: 



-4tt 2 E V 

12(1 - v 2 ) b 


(8) 


Instability of a single simply supported bead cross section is 
the final mode to be analyzed. Consideration must be made for the 
critical axial compressive stress as well as the critical compressive 
stress due to bending. The equations which define this mode (ref. 2) 
are respectively: 



Angle 02 t degrees 


Figure 29. Effects of varying angle Qz on the 
diagonal mode of instability . 
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F cb = -0.82 E(t/R) 


(9) 


1.19 

F Bb = -0.77 E{t/R) ] * 15 (10) 

On the basis of these equations, the critical buckling loads for the 
beaded panel were determined and are presented in table 7. 


Table 7. Semiclassically derived bucklmo loads 


Instabil ity 
mode 

Load 

Ibf/in (N/cm) 

Total load 

General - My Cr 

1818 

34,700 Ibf 


(3200) 

(15^,350 N) 

Diagonal - Ny DCr 

2240 

42,700 Ibf 


(3920) 

(190,000 N) 

Flat - F cf . 

- - - - 

687,000 Ibf /in 2 
(473,700 N/cm 2 ) 

Bead 1 



Axial - F Cb 


343,200 lbf/in 2 
183,400 N/cm 2 ) 

Bending - F Bb 

- - - - 

372,200 Ibf/ in 2 
(256,600 N/cm 2 ) 
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COMPARISONS OF EXPERIMENTAL RESULTS 
WITH SEMICLASSICAL ANANLSIS 

Comparisons can be made between the results of experimental test- 
ing and semiclassical analysis. As previously discussed, the critical 
elastic buckling load derived from semiclassical analysis was 1,818 
LBf/in (3,200 N/cm), the critical mode being general instability. 
However, the semiclassical analysis ignores the restraint at the ends 
of the panel which the doublers provide. Furthermore, the semiclassi- 
cal analysis assumes a constant cross-sectional thickness along the 
panel length. 'Because of these simplifications, it is probable that 
the actual buckling load should be higher. The results of the force/ 
stiffness analysis support that assertion. From figures A4(a), (b), 
(d), (e), and (f) the indicated buckling load for general instability 
was found to be between 2,660 to 2,880 LBf/in (4,660 to 5,040 N/cm) or 
about 37% higher. 

The results of the experimental analysis has shown that it was 
not general instability which was the critical mode, but a mode which 
was similar to the diagonal mode. The buckling load for diagonal in- 
stability calculated using the semiclassical analysis was 42,700 pounds 
(190,000 N). The results of the test show that the panel underwent 
elastic buckling at about 40,000 pounds (177,930 N) in a mode similar 
to that assumed by the semiclassical diagonal mode analysis. However, 
the test results show that the buckled shape of the beads did not match 
the edge conditions assumed in the semiclassical analysis (i.e., the 
edges did not remain straight). Therefore, the fairly close agreement 
of the buckling loads for this particular panel cannot indicate the 



general validity of the semiclassical diagonal mode analysis. 

The inadequacies of the semiclassical analyses, point out the 
need of more sophisticated analysis such as Nastran or more realistic 
semiclassical theory. 
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SUMMARY OF RESULTS 


For purposes of comparison, the primary experimental results will 

be reiterated at this point along with a summary of the results of the 

finite element and semi classical analysis. 

Experimental Results 

1. The panel was eccentrically loaded (due to panel geometry) as 
verified by the experimental data. 

2. Elastic buckling occurred at a load of about 40,000 pounds 
(177,930 N). The critical axial force per unit width of the 
panel was 1,680 LBf/in (2,940 N/cm). The elastic buckling load 
involved out-of -plane deformation in addition to a lateral de- 
flection characteristic of the .diagonal mode of instability. 

3. The ultimate strength of the panel/side stiffener assembly was 
48,600 pounds (216,180 N). 

Finite Element and Semi classical Analysis Results 

1. An elastic buckling load of 31,700 pounds (141,000 N) was cal- 
culated using finite element model EDGE2, 26% lower than the 
actual failure load. 

2. An elastic buckling load of 50,000 pounds (222,400 N) was cal- 
culated using finite element model EDGE3, 25% higher than the 
actual failure load. 

3. Both of the models used for static analysis exhibited an inability 
to account for the nonlinear effects of out-of -plane bending 

with loads applied in the plane of the panel. 

4. The results of the finite element buckling analysis done with 



models BEAD, FLAT and DIAG, compared poorly with the semi classi- 
cal results. 

Elastic buckling loads of 34,700 and 42,700 pounds (154,350 and 
190.000 N) were calculated for the general and diagonal modes of 
instability, respectively. These values are based upon simplified 
semi classical theory. 



CONCLUSIONS AND RECOMMENDATIONS 


The semi classi cal analysis used to design the beaded panels is 
based upon some simplifying assumptions. First, the sections used 
in the semi classical analysis were assumed to be simDly supported. 
However, the test results show that the buckled shape of the panel 
did not match the assumed edge conditions used in the semiclassical 
analysis. Secondly, significant restraint is ignored in the semi- 
classical analysis by assuming a constant cross-sectional thickness 
along the panel length. Therefore, the fairly close agreement in the 
buckling loads, for the diagonal mode in particular, cannot indicate 
the general validity of the semiclassical analysis. 

The experimental results suggest that the semiclassical analysis 
is conservative. It can therefore be concluded that the panel was 
conservatively designed and thus heavier than necessary. 

On the basis of the results of the semiclassical and finite 
element analysis the following recommendations are made: 

1. The inadequacies of the semiclassical analysis point out the need 
for more sophisticated analyses. The analysis should include 
realistic edge support assumptions as well as the use of accurate 
cross-sectional thicknesses in the analyses. 

2. The results of the finite element buckling analysis compared 
poorly with the experimental and semiclassical results. The 
reason for the poor comparisons is not precisely known, but is 
certainly an area for further research. One possible area for 
investigation would be the development of a full panel model for 
buckling analysis. 
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Finally, several concluding remarks should be made about the test 
monitoring techniques used during the tests. The force/stiffness tech- 
nique provided unconservative predictions of the elastic buckling 
strength of the panel. Furthermore, for the best results the gages 
should be in close proximity to the location of the elastic failure. 
Secondly, the Moird fringe ‘technique proved to be extremely useful as 
an aid in identifying the mode shapes of the panel. Furthermore, the 
MoirS fringe technique made it possible to identify mode shape changes 
which might have gone unnoticed based upon strain gage results alone. 
The technique also provided an accurate means of measuring out-of-plane 
displacements of the entire panel. 
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APPENDIX A 
EXPERIMENTAL DATA 

Strain gage locations are indicated on all of 
strain gage plots presented in this section, by 
two symbols. A closed circle ( • ) is used to 
indicate qages on the side of the panel shown. 
An open circle ( o) is used to indicate gages 
on the opposite side of the panel. 
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( a ) 36,000 pounds (160,140N) load. 


(d) 40,000 pounds (177,930N) load. 


Frigate M. ( continued ) Vo-te di^esienceA in cuAvaluAe 

centeA two ^latA the panel in (c) and Id). 







(e) 42,000 pounds (186,825N) load. 

Figure Al . 
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(f) 44,000 pounds (195,720N) load. 
( continued ) 






(g) 48,500 pounds (2l5,750N) load. 

Figure Al . 


(h) The panel after failure , no load. 

(continued) 
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(i) 

Figure Al . (continued) The panel after failure 
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(continued) 
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The panel after failure. 
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Figure AS. Continued 
Solid symbols denote DT data. 
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Figure A2. Continued 
Solid symbols denote DT data. 
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Figure A2. Continued 
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Figure A3. Continued 
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Figure A3. Continued 
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Figure A3. Continued 








Compressive strain, yin/in or y cm/cm 


Compressive load, lbf 


0 10 20 30 40 50xl0 3 



1 I I ! 1 I 

o 50 100 150 200 250x10 

Compressive load, n 

(q) Strain gage 426. 

Figure A3. Continued 
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Figure AS. Continued 
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Figure AS. Continued 
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Figure A3. Continued 
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The Force/Stiffness Technique 
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The Force/Stiffness Technique 


The force/stiffness technique is a nondestructive test technique 
used to experimentally determine the buckling strength of a structure 
(ref. 13). This method of nondestructive, testi no is based upon the out- 
of-plane deflection characteristics of a structure under compressive 
load. 

In figure Bl, two strain qages are shown mounted to opposite sides 
of a panel which is under a compressive load. Initially, as the load 
is increased, both qaaes measure a comoressive strain. As the column 
deflects to produce the stresses indicated, qaqe A measures an addi- 
tional compressive component due to bending and gage B measures an 
additional tension component due to bending. When the output of gage 
B is subtracted from gage A, the resulting strain is that due to Panel 
bending, only. When the compressive load is divided by the difference 
of the gages and plotted aqainst load, the result is a plot similar to 
that shown in figure Bl . Theoretically, buckling occurs when the curve 
intersects the load axis. Usually the loading is stopped before the 
curve intersects the load axis and the curve is extrapolated to an in- 
dicated (predicted) buckling load. 
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Figure B 7. Mechanics of the force, /stiffness technique . 
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APPENDIX C 

The MoirS Fringe Technique 


121 



The Moird Fringe Technique . 


The second method of measuring out-of-plane displacements of the 
panel utilized the Moird frinqe technique (ref. 14). This technique is 
best described with the aid of figures Cl and C2. The Moird frinqe 
technique requires the use of a camera, a noint source of light, and a 
grid plane, arranged in a manner similar to that shown in figure Cl. 

The grid plane referred to above is typically made on a sheet of 
clear distortion free plastic photographic film. A system of equally 
spaced parallel black lines is then applied to the film in densities 
ranging up to 500 lines/in (200 lines/cm). The greater the line den- 
sity-, the greater the sensitivity to out-of-plane displacement. 

When light is Dassed through the grid plane, shadows of the lines 
are cast upon the test specimen as shown in figure C2. As the test 
specimen deflects out-of -plane (i.e., moves either toward or away from 
the grid plane) the shadows appear to move creating fringes of dark and 
light areas. A dark fringe is formed when the shadows from the grid 
plane fill the spaces between lines on grid plane. Light areas occur 
when the shadow falls directly beneath a grid plane line. A calibra- 
tion photo of the beaded panel with no load is shown in figure C3. 
Proceeding from a dark fringe throuqh a light fringe to another dark 
frinqe in this fiqure is equivalent to an out-of-plane distance of 
0.048 inches (0.122 cm). Therefore, by selecting a stationary reference 
point (for all load conditions) the out-of -plane dimensions (displace- 
ments) can be determined. The stationary point selected in this case 
was the bottom of the Moird fringe glass assembly where the assembly 



was rigidly attached to the panel. 

The fringe in the down beads shown in figure C3 became indecern- 
able as the beads approached the flats. Therefore, paper strips., pointed 
out in figure C3 were glued to the panel as shown, in such a manner that 
the fringes could be counted down to the peak of the bead. Out-of- 
plane displacements of the down beads could then be determined. 
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Sensitivity (<$) = 1/(P Tan ©) mches/frmge 
Where P is the grid pitch in lmes/in 


Figure C2. Meehanios of the Moire fringe technique. 
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Figure C3. Sample Moire fringe calibration photograph of the beaded 

panel. 
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THE DOCTOR OF ENGINEERING PROGRAM 


The Doctor of Engineering (D.E.) program at the University of 
Kansas is a mission oriented program designed to develop creative 
leadership abilities in its participants. The D.E. program is char- 
acterized by a less specialized, broader curriculum than the typical 
Doctor of Philosophy program. Part of the D.E. curriculum consists 
of core courses including engineering management, finance, design, 
systems engineering, and sociotechnological interaction. 

The core courses are utilized by the D.E. candidate to develop, 
plan, and manage an engineering project. The D.E. proj'ects are of 
sufficient complexity to require several team workers, activity 
scheduling, and resource management. The technical and managerial 
results of the project are then reported in the form of a doctoral 
dissertation. 



PROJECT SCHEDULING 


Although this project could not be considered extremely large, 
it was of sufficient complexity to require activity scheduling and 
coordination to run smoothly. Test preparation as well as the actual 
testing of the panel required the greatest degree of coordination 
since a larger number of people were involved. 

Table D1 itemizes the major activities of this project and their 
respective durations. One item not included in Table D1 is the ac- 
tivity related to the development and utilization of the finite element. 
Nastran model, EDGE3. EDGE3 was not developed until late in the pro- 
ject (after experimental testing) and was therefore, not included in 
the original scheduling. Figure D1 is an arrow diagram indicating the 
required sequence of activities as well as the critical path. Figure 
D2 is the corresponding timetable of activities. As indicated by 
figure D2, total project duration was about 300 working days or about 
14 months. Note that neither the arrow diagram or the timetable in- 
cludes time spent at the University of Kansas attending classes full 
time. 
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Table 01: Major Project Activities 


Activity 

Remarks 

Duration 

Days 

A 

Define grid point geometry, elements, etc. 
for Nastran models EDGE1 and EDGE2. 

: 

21 

B 

Prepare case input for models EDGE! and EDGE2 
to design edge supports. 

28 

C 

Prepare case input to model EDGE2 to analyze 
general instability mode. 

28 

D 

Define grid point geometry, elements, etc. for 
Nastran models BEAD, FLAT, and DIAG. 

14 

E 

Analyze other instabilities with BEAD, FLAT, 
and DIAG with appropriate input. 

28 

F 

Design and fabricate panel end supports and 
side stiffeners. 

42 

G 

Design and fabricate Moirg grid glass supports 

7 

H 

Draw up data aquisition system plans and set 
up data aquisition system. 

21 

I 

Assemble edge supports to panel. 

2 

J 

Machine assembled panel ends parallel. 

3 

K 

Draw up instrumentation plans and instrument 
panel 

14 

L 

Set up testing machine, photographic fixtures 
and other equipment. 

7 

M 

Adhere Moirg grid to glass and assemble glass 
supports and glass. 

7 

N 

Establish loading rates and make preliminary 
check of system. 

1 

0 

Test panel. 

21 

P 

Analyze data. 

28 

Q 

Prepare report. 

150 
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Figiive D2. Timetable of maj 
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RESOURCE MANAGEMENT 


The majority of this project was completed at NASA's Dryden 

* 

Flight Research Center. As a result, a number of NASA personnel 
were involved with this project to some degree. The NASA personnel 
were directed in their activities (for this project), by the princi- 
pal investigator (naturally through the appropriate channels). 

Table D2 is a list of the personnel involved with this project and 
the time required of each. 

The greatest contribution, both in terms of time and consultation, 
came from Roger Fields, an aerospace engineer in the aerostructures 
division of the research directorate at DFRC. Mr. Fields was the 
principal investigator's immediate supervisor at DFRC and an ad hoc 


Table D2: Manpower Requirements 

Person 

Total Hours 

Principal Investigator 

2000 

Ad Hoc Committee Member 

750 

Committee Chairman 

200 

Operations Engineer 

20 

Electronics Technicians- (2) 

120 

Instrumentation Technician 

50 

Photo Lab Personnel 

20 

Technicians- (2) 

30 

Student Aids 

40 

Machimsts-(3) 

30 

Total Man Hours 

3260 



member of the committee supervising this project. Mr. Fields spent 
countless hours reviewing design procedures, experimental data, and 
ultimately this report. His advise was particularly useful from a 
practical (applications) standpoint. 

Other NASA personnel involved with this project include Gary 
Craton and Bill Hensley, electronics technicians responsible for the 
setup and operation of the data acquisition computer and related 
electronics. Walter Sefic and Larry Reardon, operations engineers, 
were responsible for the operation of the MTS testing machine during 
each test. Clarence Cook, an instrumentation technician, instru- 
mented the panel with strain gages. Two operations technicians, 

Ralph Sparks and Leon Hatcher, prepared the MTS machine for use and 
installed brackets on to which displacement transducer and photo 
equipment were mounted. Other personnel included machinists, ohoto 
lab technicians and student aids who reduced and analyzed data. 

Technical and theoretical direction was provided by Dr. John 
Easley, professor of civil engineering at the University of Kansas. 

Dr. Easley, spent long hours reviewing the paper for correctness of 
its assertions and conclusions, as well as making two trips to DFRC 
to observe testing procedures and make suggestions. Dr. Easley was 
also the chairman of the committee supervising this project. 

Another asset to this project was the availability of NASA’s 
physical facilities, test equipment, and computer facilities. Many 
hours of finite element computer programs were run on the CDC computer 
at DFRC. Had the CDC computer not been available, it is doubtful that 
the computer analysis done for this project could have been accomplished. 



The Nastran model EDGE3 was so large, however, that it was necessary 

to have the buckling analysis done with EDGE3 run on the University 

of California at Los Angeles' IBM computer. Other computer facilities 

included a SEL 810 data acquisition and control computer complete 

with peripheral devices such as TV -momters, a line printer, and strip 

charts. The testing machine used to test the panel was an MTS machine 

capable of programable loading rates. Other facilities previously 
* 

mentioned include an instrumentation lab, a machine shop, and a pho- 
tography lab. 
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